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Abstract 



Srivastava et al. ( 2010() have observed a highly twisted coronal loop, which was 



anchored in AR10960 during the period 04:43 UT-04:52 UT on 4 June 2007. The 
loop length and radius are approximately 80 Mm and 4 Mm, with a twist of 
11.5 7T. These observations are used as initial conditions in a three dimensional 
nonlinear magnctohydrodynamic simulation with parallel thermal conduction 
included. The initial unstable equilibrium evolves into the kink instability, from 
which synthetic obscrvables are generated for various high-tcmpcrature filters of 
SDO/AIA. These observables include temporal and spatial averaging to account 
for the resolution and exposure times of SDO/AIA images. Using the simulation 
results, we describe the implications of coronal kink instability as observables 
in SDO/AIA filters. 
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1. Introduction 



The active region magnetic field exhibits complex topology due to the 
emergengce of new large-scale magnetic flux from the sub-photosphcric 
layers of the Sun. Such complexity some times generates various types of mag- 
netic instabilities that are well observed in the solar active regions as a trigger of 
flares and associated dynamical processes ( [S rivastava et al. . 2010l ; Kumar et al 



2010t Foullon et al. , 2011 ; Innes et"all 120121) . Among the various types of mag- 



netic instabilities, the kink instability is most commonly observed in the solar 
atmosphere, and evolved in the large-scale magnetic flux-tubes due to their az- 
imuthal twist. However, the twist must cross the minimum threshold value of 
$ > 2.5 7T to trigger the kink instability in the magnetic flux-tubes in the solar 



corona (IHood fc Priestl. 11979 ) . The kink instability is well observed in active 
region loops ( Srivastava et al. . 2010)) . eruptive filaments ( Liu et al. ■ l2008h . erup- 
tive cor onal cavities (Liu et all 120071). an d it ca n eventuall y trigger flares and 



CMEs (jKliem fc Torok 1 . 120061: Icho et all [20091) . Recently. [Zaaarashvili et al 



( 20101 ) have given an important conclusion that the axial mass flow in the mag- 
netic flux-tubes reduces the threshold of kink instability, and can easily lead to 
solar flares and CMEs. The kink instability does not only liberate the flare en - 
ergy, it can also cause the excitation of MHD wave modes (jHavnes et al. . 2008). 
The kink instabilities observed in the as trophysical plasma are also firstly re- 
produced in the laboratory experiment (jMoser fc Bellanl . I2012I ). which shows 
the confinement of kink unstable twisted plasma structures due to multi-stage 
cascade reconnection at the laboratory scales. 

In addition to the observations the coronal kink instability is studied using 
analytical and numerical modelling dMikic et al. , 1990t iBatv fc Hevvaertsl . 1996 



analytical an d numerical modelling (Mikic et aUIHJiJUtl-oatv &; ri evvaerts, lltjbifr ; 
Arber et all 1 19991 : iGerrard et ai] . l2002l ). The ideal MHD kin k instability has 



been inferred as a trigg er for reconnection in the flaring loops ([Browning et all 
l2008HHood et all 120091 ). iBotha et all ([201 1[ ) initialised the kink instability with 
a twisted magnetic field and reported that the inclusion of thermal conduction 
along magnetic field lines can reduce the maximum reconnection generated tem- 
perature by an order of magnitude. This changes the range of spectral lines as 
observab les of the kink inst ability when compared to simulations without con- 
duction. IBotha et all (|2012h have also reported the observational signatures of 
the MHD kink instability with the inclusion of thermal conduction. A kink- 
unstable cylindrical loop is evolved b y using the initia l cond itions of a highly 



twisted loop as observed by TRACE ( Srivastava et al. . 2010l ). The numerical 



results arc synthesized through TRACE temperature response function filtering, 
and through spatial and temporal averaged line of sight intensity measurements 
( Havnes fc Arberl . 120071 ). 

In the present paper, we extend the work of IBotha et al. I d2012l) and stu dy 
the evolution of the kink unstable loop as observed by Srivastava et al. ( 2010l ) in 
the high-temperature filters of SDO / AIA. In Sec. 2 we present a brief description 
of the numerical model and initial conditions. In Sec. 3 we describe the results 
related to the synthesis of a kink unstable loop in various SDO/ AIA channels. 
The discussion and conclusions form the last section. 
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2. Numerical Simulation of the Observed Coronal Kink Instability 



The nonlinear three-dimensional simula tion is executed by implementing the 
MHD Lagrangian-remap code (Lare3d) of lArber et al.l (|200lh . The numerical 
code solves the resistive MHD equations for the fully ionised plas ma with a 



heat flux included in the energy equation (jBotha et al.l . l201ll 120121 ) . Thermal 



conduction is considered to be incl uded along the magnetic field lines in the form 
of a classical ISpitzer-Harml (|1953l ) or Braginskii conductivity with log A = 18.4 
that corresponds to the standard therm al conductivi t y parallel to t he ma gnetic 
field of «|| = lCT 11 ! 15 / 2 W m" 1 KT 1 (|Priestl 1200(1 Botha et ail . l2012t ). The 
numerical code consists of an artificial resistivity that is activated only when 
the electri c current exc e eds a critical value, which is set in the present case to 
j c = 2 mA (jBotha et al. 



as follows 



20121 ). This resistivity is considered in the form given 

(1) 



o. 



\j\ > Jc 

lil < jc 



where ryo is the anomalous resistivity. rjQ is switched on when the kink insta- 
bility occurs and thereaft er it remains active for the duration of the simulation 



( Botha et al.l . 1201 ll 120121) . The coronal loop is initialised as a straight twisted 
cylinder in a uniform background temperature and density (Fig. 1) as previously 



performed bv lBotha et~alT(|2012l) . 

The initial conditions of the numerical si mulation of the coronal l oop is taken 
from observations of a highly twisted loop ( Srivastava et al.l . l2010l ) . Using the 
multi- wavelength observations of SOHO/MDI, SOT-Hinode/blue-continuum (4504 
A), G band (4305 A), Ca II H (3968 A), and TRACE 171 A, they observed a 
highly twisted magnetic loop in AR 10960 during the period 04:43 UT-04:52 UT 
on 2007 June 4. SOT-Hinode/blue-continuum (4504 A) observations show that 
the penumbral filaments of a positive polarity sunspot have counterclockwise 
twist that may be caused by the clockwise rotation of the spot umbrae and can 
activate a right-handed helical twist in the loop system. This loop whose one 
footpoint is anchored in this sunspot, shows strong right-handed twist in chro- 
mospheric SOT-Hinode/Ca II H (3968 A) and coronal TRACE 171 A images, 
which was consistent with the Hemisphere Helicity Rule (HHR). The length and 
the radius of the loop are estimated as L«80 Mm and a«4.0 Mm, respectively. 
The distance between neighbouring turns of magnetic field lines is estimated as 
wlO Mm, which further gives the total twist angle $=11.57r as estimated for 
a homogeneous distribution of the twist along the loop. This observed twist 
is much larger than the Kruskal-Shafranov instability criterion (<I>>27r) for the 
kink instability and later triggers a B5.0 class solar flare that occurred between 
04:40 UT and 04:51 UT in this activ e regio n. The details of these observational 
finding are given in ISrivastava et al. I (|2010l) . 
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Figure 1: Left-panel: The cartoon shows an orientation of the cylindrical loop in the Cartesian 
numerical domain. The dotted line indicates about the integration path of the line of sight 
integral as described by Eq.3 along the y-direction and perpendicular to the x — z plane. 
Right-panel: Temperature and mass density profiles along the length of the cylindrical loop 
on its central axis at a time of 579.98 s. The solid (dashed) line is the temperature (density). 
Initialisation is occurred with a constant temperature of 0.125 MK and a constant density of 
1.67 ng m~ 3 . 



The observed curved loop is straightened to perform numerical 
simulations (Figure 1, left-panel). One of t he footpoints of the a ctive 
region loop is anchored in a positive polarity spot ( Srivastava et al. . 2010t ) and 
it is es timated that the m inimum value of the magnetic field at this location is 
470 G (jBotha et all 120121 ) . From this photospheric value of the magnetic field, 
the chromospheric field is calculated approximately as 80 G, as well as near 
the l oop apex in the corona as 20 G (|Petrie fc Patrikeeval [20091: iBotha et all 
20121) . The numerical simulation is initialised for the loop with a maximum 



internal field of 20 G, while the outside background magnetic field is considered 
as a uniform value of 15 G parallel to the cylindrical axis of the simulated 
loop structure. In the numerical simulations the coronal loop is initialised as 
a uniform cylin der in force- free equilibrium and set u nstable to an ideal MHD 
kink instability (jHood et all l2009t IBotha et all 120121) . The axial twist is given 
by the following equation 



LB e 

rB z 



with max($) = 11.57T, 



(2) 



where L is the loop length and r is the radial distance from the central axis of the 
loop. The axial and azimuthal magnetic fields are given respectively by B z and 
Bg as a function of r. The maximum twist is considered at position r = 1 Mm, 
which exceeds the stability threshold , therefore, the loop becomes kink unstable 
dMikic et al. . 1990f ) . The details of the magnetic field configuration is outlined in 
Botha et all (|201ll . l2012l) . Gravity is not included in the simulations. The initial 



temperature and mass density arc taken as uniform and constant, with a typical 
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time [s] 



time [s] 



Figure 2: Left-Panel: Maximum temperature during the evolution of the kink instability 
in which the solid (broken) line is without (with) thermal conductivity. Right-Panel : The 
average temperature during the evolution of the kink instability in which the solid (broken) 
line is without (with) thermal conductivity. The minimum temperature is also represented by 
the dash-dot-dashed (dotted) line without (with) thermal conduction. 



values of 1.67 x 1CT 12 kg m" 3 (jYoung et all 120091 ) and 0.125 MK respectively. 
This temperature is selected to make the evolution of the kink instability as 
evident in the SDO/AIA filters. During the evolution of the kink instability the 
temperature i ncrea ses locally where reconnection occurs in the simulated loop 
( Botha et all l20lH) . This temperature increament is sufficient that the results 
are highly insensitive to the lower initial temperatures. The initial temperature 
was chosen specifically so that these high temperatures (Figs.[TJ right-panel and 
\2§ will be visible using the SDO/AIA filter response functions (Fig. [3]). The 
details of these tec hnical aspects of the numerical simulation are described in 
iBotha et all (j2012t) . 



3. Synthesized Kink Unstable Loop as an Observable in Various Fil- 
ters of SDO/AIA 

The evolution of the kink unstable coronal loop is studied in the Cartesian 
geometry (Fig. 1) of the simulation domain. The boundaries in the x — y plane 
perpendicular to the loop axis are considered at ±8 Mm and it is set to be 
reflective. Therefore, there are 4 Mm seperation between the loop edge and 
the outer boundaries to avoid the numerical reflections near the boundaries. 
Along the loop axis, these boundaries are cosidered at ±40 Mm with velocities 
and temperature fixed respectively at zero and the initial background value by 
allowing the temperature gradients. Therefore, the heat flux across the ends of 
the loops exists in the simulation domain. The grid resolution in x — y — z is 
given by 128 x 128 x 256. These technical details are summarized in lBotha et al.l 



(12012ft 



The images generated from the numerical simulations were obtained by using 



the modified temperature response functions of SDO / AIA (jAschwanden fc Boerner 
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2011). Details of the line contributions for the AIA channels are outlined in 



O'Dwver et al. also. The emission is calculated at every node of the 

numerical grid and thereafter integrated along the y-direction perpendicular to 
the axis of the simulated loop (cf., Fig.[TJ left-panel). The line of sight intensity 
integral is given by 

V 9(T)p 2 dy (3) 



where /, g(T), p are respectively the measured intensity, the temperature re- 
sponse function of the various AIA filters, and the mass density. Fig. [3] presents 
the temperature response functions of all AIA filters. This integration produces 
intensity images in the x — z plane, which are then integrated over time with 
the typical time interval determined by the exposure time for SDO/AIA as 2.9 
s. After this procedure, the time integrated image is degraded by spatially aver- 
aging over squares of 0.375 x 0.375 Mm 2 to compensate for the pixel resolution 
of SDO/AIA. The similar process is adopted for the generation of synthesised 
TRACE observables as reported bv lBotha et al. ( 2012 ). However, in the present 
study, we perform this image synthesis of the kink unstable loop as an observ- 
able in various filters of SDO/AIA (cf., Figs. 4-6). It should be noted that 
we used corrected AIA 94 response function with correction for en- 
hanced response at log(T)<6.2 due to missing Fe IX, Fe XII lines in 
CHIANTI. Similarly, the response function for AIA 131 A is used 
with considering the contribution of Fe VIII and Fe XI. This may 
vary overall intensity features in the synthesized images of 094 A and 
131 A compared to those generated by using old response functions. 
The 94 A , 211 A and 131 A images are synthesized by choosing the 
elimination of lowest 1% of the intensity values at time 321.9 s. 

After initialisation, the simulated kink unstable loop evolves through a linear 
phase with a duration of 300 s (Fig. [5]). During the nonlinear phase, the kink 
instability drives magnetic field into a current sheet where reconnection occurs 
and the temperature reaches to a maximum value due to the released energy 
and subsequent heating. Thereafter, thermal conduction along the magnetic 
field lines transports the generated heat flux along the field lines away from the 
reconnection sites where the plasma was heated. The time-scale for the evo- 
lution of the kink instability in all SDO/AIA filters is the same as of the case 
without thermal conduction. The physical processes during the kink instability, 
wi th and witho u t ther mal conduction, are disc ussed in comprehensive manner 
bv lBotha et all (fioTTh and lBotha et akl <20l3 ). In the present numerical sim- 
ulations, the thermal conductivity causes the spread of temperature along the 
magnetic field lines resulting in the synthesized images in which the features are 
less resolved. 

The synthesized images show that the footpoints of the loops increase their 
emission during the nonlinear phase of the kink instability irrespective of the in- 
clusion or exclusion of thermal conductivity. The line of sight intensity integral 
(Eq.3) is determined by the temperature response function of SDO/AIA and 
the chosen mass density. Fig. [T] (right-panel) gives the temperature and mass 
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Figure 3: Normalised SDO/AIA temperature response functions with the minimum and maxi- 
mum temperatures during the simulation represented by vertical dotted lines. The 304 A band 
captures chromospheric and transition region temperatures, features that are not included in 
our model. Consequently it is not included with the other filters. 
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Figure 4: AIA 94 A simulations with exposure time 2.9 s and spatial resolution 0.375 Mm 
per pixel. The intensity is in reverse colour, with white the lowest value on the scale. This 
minimum is chosen so that the lowest 1% of the values at time 321.9 s are eliminated from 
these plots. All the exposures use the same reverse colour scale. 



density profiles for the numerical simulation including thermal conduction. It 
is noticiable that the enhanced emission is due to a density increase at the foot- 
points of the simulated loop and it is not due to the footpoint heating. Plasma 
flows from the middle of the numerical domain where the current sheets are 
formed as reconnection sites , driven by MHD ponderomotive forces generated 
during the evolution of the kink instability. Footpoint brighte ning due to com- 
pressi on w ere also evident in the coronal loop simulations of lHavnes fc Arber 
( 20071 ) and lBotha et al.l (|2012f) . The factor that determines the initiation of the 
nonlinear phase of the kink instability is the initial twist (Eq.2) in the coronal 
loop. A smaller value of maximum twist increases the duration of the linear 
phase and vice versa. However, it is found that once the nonlinear phase is 
evolved, the formation of the current sheet, the reconnection, and the thermal 
aftermath have t he same duration as long as the initial twist exceeds the sta- 
bility threshold (IBotha et al. . 2012). Th ese physical and technical details are 



discussed in detail by IBotha et al.l ( 20121) . 

During the study of the evolution of such kink instability and its comparison 
with the observations, it is more realistic to compare the obtained results when 
thermal conduction is included in the numerical simulation model. The first 
snapshot in Fig. 2]is synthesized at 321.9 s for AIA 94 A and shows the twisted 
structure of the current sheet formed at the initiation of the nonlinear phase of 
the evolution of kink instability. In the later images, the thermal conduction has 
transported the reconnection generated heat flux from the internal reconnection 
sites along the magnetic field lines. As the internal reconnection occurs within 
the simu lated loop, the magnetic field lines become straight along the length of 



the loop (jHavnes fc Arbeit 120071 ; IBotha et al.l . l2012l ). At the same time, the heat 



flux is transported along magnetic field lines away from the reconnection sites. 
Therefore, the loop becomes cool and its internal magneto-plasma structures 



are less resolved ( Botha et all 12012) 



The corrected and modified temperature response functions for the various 
wavebands (filters) of SDO/AIA arc presented in Fig. [3] The emissions in 171 
A for SDO/AIA and TRACE are almost the same. The only difference between 
these two observational instruments is the time resolution. The TRACE ex- 
posure time is approximately eleven times longer than th at of the SDO/AIA . 



Therefore, the images from SDO/AIA are much sharper. iBotha et al.l (|201 2| ) 
have already reported the kink unstable loop as an observable in TRACE 171 
A filter, therefore, we do not present the similar results as obtained in AIA 171 
A filter. In the present paper, we report the evolution of observed kink unstable 
loop as an observable in the high-temperature filters of SDO/AIA. Moreover, 
AIA 193 A and 211 A filter responses are almost similar, therefore, we present 
only the results related to the 211 A synthesized AIA images. The 304 A band 
captures chromospheric and transition region temperatures and related features 
that are not included in our model. 

The emission images for SDO/AIA 94 A, 131 A and 211 A are presented 
in Figs. HUni As each temperature response function (Fig. [3]) is different from 
the others, the structures visible by each waveband are slightly different. The 
temperature response function for SDO/AIA 94 A captures the top half of 
the temperatures obtained during the simulation, so that the heated magnetic 
structures inside the loop are visible in the intensity images, together with the 
average temperature inside the loop (Fig. 2]) . The temperature response func- 
tion of SDO/AIA 131 A captures the middle of the temperature range during 
the simulation, so that the intensity images (Fig. [5]) shows mainly loop structure 
at the average temperature. The images obtained for SDO/AIA 211 A (Fig. [6]) 
captures the maximum temperature during the simulation, with a contribu- 
tion from the average and lower temperatures, as shown by its temperature 
response function (Fig. [3]). The emission images of SDO/AIA 304 A (Fig. [3j 
are not presented in this study. This waveband captures the temperatures of 
the chromosphere and transition region, which are not included in our coronal 
loop model. Also, the emission images obtained for SDO/AIA 335 A are not 
presented here. This temperature response function is broad and spans all the 
temperatures generated in the simulation. As a result, 70% or more of the 
lower values of the intensity obtained by (Eq.3) have to be eliminated in order 
to produce intensity plots in which the loop structure is clearly visible. 

4. Discussions and Conclusion 

A highly twisted and kink unstable coronal loop is studied numerically by 
solving the resistive MHD equations for a fully ionised plasma and with the in- 
clusion of the thermal conduction. The loop is initialised as a straight cylinder 
with a very high twist above the stability threshold leading to the kink insta- 
bility. The simulations were initialised with phys ical parameters derived from 
the observations of a highly twisted coronal loop ( Srivastava et al. . 2010l ). The 



initial magnetic field structure used in the simulations is consider ed as a kink 



unstable force-free equilibrium where the twist varies with radius (jHood et al 




, [Km] » M , [Mm] « M ■ [»">] 



Figure 5: AIA 131 A simulations with exposure time 2.9 s and spatial resolution 0.375 Mm 
per pixel. The intensity is in reverse colour, with white the lowest value on the scale. This 
minimum is chosen so that the lowest 1% of the values at time 321.9 s are eliminated from 
these plots. All the exposures use the same reverse colour scale. 



20091 : iBotha et all 120121 ) Line of sight emission intensities were calculated 



from the simulation data using the corrected temperature response functions 
from SDO/AIA. The different wavebands from the SDO/AIA exhibits different 
magnetic field structures, as each filter temperature response function captures a 
different plasma temperature and related emissions. Therefore, it becomes plau- 
sible to study many high-resolution features of a highly twisted kink-unstable 
coronal loop as observables of the various filters of SDO /AIA. During the evo- 
lution of the kink instability, the internal structure of the synthesized loop is 
shown in the generated images (Figs. 0]|6]). This structure is transformed into 
a simpler magnetic configuration as the kink instability triggers multiple recon- 
ncctions that cause the straightening of the internal magnetic field. The kink 
instability heats the plasma near the current sheet due to the reconnection. This 
heat flux is transported along the magnetic field lines away from the reconnec- 
tion sites due to thermal conduction. Therefore, the temperature maximum 
in the simulation is lowe r compared to the case when thermal conduction is 



excluded from the model (IBotha et all 120121) 



The observations by Srivastava et al.l "( 20101 ) clearly shows that after the 
heating as well as activation of high twist ,the loop plasma is condensed down 
within few minutes and becomes invisible in coronal filters. The duration of 
the simulations exhibits 4.25 minutes of the nonlinear evolutionary phase of 
the kink instability. During this time span the average temperature of the loop 
increases (Fig. [2j right-panel) while the maximum temperature with the thermal 
conduction is almost unchanged (Fig. [2j left-panel). The reason may be that 
the current sheets evolve through the multiple reconnection process instead 
of continuous h eating episode. T herefore, the heating location moves in the 



simulated loop ( Hood et al. . 20091) . Footpoint brightening due to compression 
of the plasma is also evident in the simulation , however, it is not evident from 
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Figure 6: AIA 211 A with parameters and times as in Figure[4] The minimum (white) of the 
intensity scale is chosen so that the lowest 1% of the values at time 321.9 s are eliminated 
from these plots. All the exposures use the same reverse colour scale. 



the co ronal observations by Srivastava et al. ( 2010l ). Although, Srivastava et al 



( 201dh have observed brightpoints at chromospheric temperatures that may be 
the evidence of localized footpoint heati ng in lower a t mosp here. 

In conclusion, we extend the works of lBotha et al.1 (|2012l ) and show that how 
the kink unstable observed loop will be synthesized in various high-temperature 
filters of SDO/AIA. The physical scenario and the technical details of the nu- 
merical simulation presented in the paper are very similar to the Botha ct al. 
(2012). However, the present work demonstrates about the synthesis of the ob- 
servational signature of coronal kink instability in the various filters of SDO /AIA 
and discuss its physical significance. More studies should be performed using 
forthcoming observations from SDO/AIA regarding the helical kink unstable 
loops that may shed more light on the stringent simulation initial conditions, 
and can provide more understanding of the energy generation and transport in 
the coronal loop systems during the evolution of such instability. 
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